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Abstract. Measurements of Ultra-High Energy Cosmic Rays achieved remarkable progress 
in the last 10 years. Physicists, gathered from around the world in the symposium UHECR- 
2012 held at CERN on February 13-16 2012, reported their most up-to-date observations, 
discussed the meaning of their findings, and identified remaining problems and future 
challenges in this field. This paper is a part of the symposium proceedings on the experi- 
mental summary and future prospects of the UHECR study. 

1 Introduction 

Cosmic Rays (CRs) are high energy radiations, mostly made up of protons and nuclei, arriving at the 
Earth isotropically from outer space. Their energies range from 10^ eV to 10^' eV, with a steeply and 
monotonously decreasing flux approximately proportional to E"^. Ultra-High Energy Cosmic Rays 
(UHECRs) reaching an energy of lO^*' eV start interacting with the Cosmic Microwave Background 
(CMB) during their propagation in extra-galactic space, produce nuclear A resonance and lose part of 
their energy in the decay of the A. In 1966, just after the discovery of the CMB, Greisen, Zatsepin and 
Kuzmine (GZK) predicted that this interaction causes a strong suppression of the CR flux at around 
1 0^*' e V (GZK cutoff) [ 1 1 . The search for this "end-point" in the CR spectrum has been one of the most 
important subjects of cosmic ray physics. 

Cosmic Rays were discovered by V. Hess 100 years ago in 1912 using an electroscope on board 
a balloon reaching 5,400 m altitude. Extensive Air Showers (EAS), particle cascades caused by the 
high energy primary CRs interacting with an atmospheric nucleus, were identified by R Auger in 1938 
using a coincidence of Geiger counters separated by up to 150 m on the ground 121 . The existence of 
UHECRs reaching 10^" eV was first reported by J. Linsley in 1962, using a network of scintillation 
counters deployed over an area of 8km2 t3J. 



2 Ultra-High Energy Cosmic Rays 

UHECRs are extremely rare; they are detected as Extensive Air Showers reaching the ground over 100 
km^ area approximately once per year. A traditional way of detecting UHECRs is by using an array 
of Surface Detectors (SDs) sparsely deployed on the ground. The energy and the aiiival direction 
of the primary CR are reconstructed from the number of charged particles and their arrival times 
measured by the SDs. One of the largest arrays using this technique, AGASA with 100 km^ covered 
area using optical fiber network, reported a possible extension of the UHECR energy spectrum beyond 
the predicted GZK cutoff' in 1998 [4 |. 

A different method of detecting EAS was discussed by A. Chudakov, M. Oda and K. Suga in 
1958-1962; they proposed to detect scintillation (fluorescence) light produced by the passage of EAS 
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particles in the atmosphere by using an imaging telescope [S]. The first detection of EAS by air flu- 
orescence was made in Japan by G. Tanahashi in 1968 |6|, who came back from Cornell University 
where an intensive study of air fluorescence detection had been initiated by K. Greisen. The method 
was systematically developed and refined by the Fly's Eye group in Utah, and the first observation 
of the GZK cutoff was announced in 2008 |7| by the upgraded detector. High Resolution Fly's Eye 
(HiRes). 

The history of EAS observation and fluorescence method was reviewed by P. Sokolsky^ at the be- 
ginning of UHECR-2012 (for a historical review of EAS, see also [8J by K. H. Kampert and A. A. Wat- 
son). Recent efforts of standardizing the air fluorescence yield and its dependence on the atmospheric 
condition were reported by B. Keilhauer. 

The emerging air fluorescence method enabled: 

- experimental determination of the primary CR energy by measuring the total amount of fluo- 
rescence light in the atmosphere (total absorption calorimetry), and 

- statistical determination of the primary particle species, or the mass composition, by measuring 
the longitudinal shower development in the atmosphere (Kmux measurement, see below). 

Modern large-scale UHECR experiments, Pierre Auger Observatory (hereafter called Auger) and 
Telescope Array experiment (TA), are hybrid experiments employing a large SD array and a battery of 
Fluorescence Detectors (FDs) overlooking the SD area. The SD array has a high duty factor (~100%) 
and has flat acceptance above a certain threshold energy. The FD is operated only on clear, moonless 
nights (duty < 10%), and its acceptance grows with the energy of CRs. The number of SD/FD coin- 
cident events (hybrid events) are limited, but they are important as they carry the largest amount of 
information on the produced EAS. For Auger and TA, high statistics analyses of energy spectrum and 
anisotropy are performed using the SD data, with its energy calibrated by the FD data using hybrid 
events. The mass composition is determined by using average X„,aj;, < X,„ax >, and its distribution 
width, RMS(X„,fl v), observed by the FD, where yi„u,x [g/cm^] is defined as the slant atmospheric depth 
of the maximum shower development in the atmosphere. 

The SD array of Auger is composed of 1600 water Cherenkov counters deployed in a triangular 
grid of 1.5 km spacing, and covers a ground area of 3000 km-. Each tank contains water in a bag 
with 10 m^ bottom area and 1.2 m height. The SD of TA is composed of 507 scintillation counters 
deployed in a mesh with 1.2 km spacing, and covers an area of 700 km^. Each TA/SD has two layers 
of plastic scintillators, 1 .2 cm thick and 3 m^ large. Water Cherenkov counters and plastic scintillators 
are sensitive to both the Electro-Magnetic (EM) component (electrons and gammas) and hadronic 
component (muons) in the EAS. The Cherenkov signal in the water tank is larger for penetrating hard 
muons than soft electrons and gammas in the shower, thus making the Auger/SD more sensitive to the 
hadronic component in the EAS. The plastic scintillator of TA equally samples the charged particles, 
and the resultant signal from the EAS is dominated by the outnumbered electrons and photons. 

The number of observed UHECRs increased by more than an order of magnitude by the new exper- 
iments. Auger and TA, and the quality of UHECR observation improved largely by the simultaneous 
use of the ground array (SD) and air fluorescence telescope (FD). A combination of the two experi- 
ments. Auger in Argentina (35.3° South) and TA in Utah, USA (39.4° North) now covers the entire 
sky for CR observations. 



3 Symposium UHECR-2012 

The symposium UHECR-2012 was held at CERN on February 13-16 20 120 as a follow up of UHECR- 
2010 held in Nagoya in December 20100 

There were 230 participants from 25 countries, and 53 oral and 48 poster presentations were made. 
Five working groups (WGs) were formed before the symposium from members of Auger, HiRes, 

' The name in italic is meant his/her contribution at the UHECR-2012 is available in the printed proceedings 
or in the form of presentation file at the symposium website. 
^ http://indico.cem.ch/conferenceDisplay.py?confId=152124 

^ http://uhecr2010.icrr.u-tokyo.ac.jp/ and http://proceedings.aip.Org/resource/2/apcpcs/ 1 367/ 1 ? | 
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Table 1. UHECR-2012 Working Groups. 



Subject 


Convenors 


Members 


Spectrum 


B. R. Dawson', 
Y. Tsunesada' 


T. Abu-Zayyad2•^ D. Ikeda\ D. Ivanov^'^, I. C. Maris', 
M. Pravdin'', M. Roth', A. Sabourov'*, F. Salamida' 


Composition 


J. Bellido', 
J. W. Belz^'^ 


E. Barcikowski', Y. Egorov"*, S. Knurenko"*, V. de Souza', 
Y. Tameda', Y. Tsunesada^, M. Unger' 


Anisotropy 


R Sommers', 
P. Tinyakov^^ , 


O. Deligny', A. Ivanov'', J. de Mello Neto', 
H. Sagawa', L. Timofeev'', I. Tkachev' 


Multi-Messenger 


M. Risse', 

G. I. Rubtsov^ * 


J. Alvarez-Muniz', B. T. Stokes-^ 


Shower Simulation and 
Hadronic Interactions 


R. Engel' 


J. Allen', A. Castellina', K. Itakura*, K. Kasahara'-^ 

S. Knurenko"*, S. Ostapchenko', T. Pierog', A. Sabourov'*, 

T. Sako^ B. Stokes2-\ R. Ulrich' 



1: Auger, 2: HiRes, 3: TA, 4: Yakutsk, 5: LHCf 
*: KEK Theory Center 



TA, Yakutsk air shower array (Yakutsk) and the very forward LHC experiment (LHCf). The physics 
subjects of the WGs are listed in Table-[T] together with WG members and convenors. 

The WGs addressed such questions as: how much common understanding have we reached in 
UHECR observations, what are the remaining differences among experiments, and what are the charges 
and challenges to resolve the issues. Each WG reported the results of their research in the symposium, 
and reports are included in the proceedings ||9] ifTOl llTTI lfT2l ifTH . Two summary talks were given at 
the end: one for the theory and phenomenology by A. Olinto, and another for the experiment and future 
prospects by M. Fukushima (this report). The symposium concluded with remarks by A. A. Watson and 
a general discussion among all participants. 



4 Energy Spectrum 

The spectrum WG reviewed the energy spectrum of cuiTently operating experiments. Auger, TA and 
Yakutsk, and of completed experiments, AGASA and HiRes. First, all spectra were fitted with broken 
power laws. The strong suppression of the flux, expected from the GZK process, was confirmed with 
good statistical significance by HiRes, TA and Auger around lo'^-^-''-^ eV, showing a decrease of the 
power law index from -2.7 down to -5.2 ~ -4.2. This suppression was not seen in the AGASA 
spectrum, and Yakutsk currently has insufficient exposure in this energy region. A dip structure, or 
"ankle", was observed by all experiments as a change of power law index from -3.3 to -2.7 around 
lO'^-' eV. 

Second, the WG examined the consistency of spectral shapes between experiments by adjusting 
the energy scale of each expeiiment to a common scale, which was taken halfway between Auger and 
TA. This particular choice of common energy scale, a mean of Auger and TA, reflects our belief that 
the energy determination by FD is experimentally better controlled than SD without relying too much 
on the result of EAS simulations and their hadronic interaction models at UHE. Fig. [T|was obtained by 
fitting the spectra below 10'^^ eV to the common energy scale. Corresponding scale factors are listed 
in Table-El |9 1 . 

The estimated energy scale uncertainties of FD measurement are 17%, 22% and 21% for HiRes, 
Auger and TA respectively, and the obtained scale factors are within these limits. The AGASA scale 
factor (0.652) is at the limit of consistency with the quoted systematic uncertainties of AGASA (20%) 
and FD (~20% typical). AGASA's number of events above 10^" eV becomes two after the energy 
rescaling, and the observation of the extended spectrum beyond the GZK cutofi" becomes statisti- 
cally insignificant. The Yakutsk SD spectrum is calibrated by the Cherenkov light emission and its 



EPJ Web of Conferences 



E[eV] 



> 



M) 
O 



24.8 
24.6 
24.4 
24.2 

24 
23.8 
23.6 
23.4 
23.2 

23 
22.8 
22.6 



2x10 



10' 



2x10' 



Auger (ICRC 2011) (Ex 1.102) 
Telescope Array (E x 0.906 ) 
Yakutsk (Ex 0.561 ) 
HiResl (Ex 0.911) 
HiRes 11 (E X 0.903 ) 



_L 



_L 



18.5 



1' 



19.5 



20 



20.5 



log^^(E/eV) 



Fig. 1. Re-scaled energy spectra in the form of x flux. See spectrum WG's report for details (9]. 
Table 2. Energy scale factors obtained by fitting the spectrum ||9)- 



Auger TA AGASA Yalcutslc HiRes-1 HiRes-lI 
1.102 0.906 0.652 0.561 0.911 0.903 



rescale factor (0.561) seems to be also at the limit of consistency with the systematic uncertainty of 
the Cherenkov calibration (~25%). 



5 Models 

Observed energy spectra are interpreted by models of CR generation and propagation. In those models, 
the generation of CRs at the source is usually characteiized by a power law spectrum and the maxi- 
mum energy of acceleration. The propagation of CR causes energy losses by collisions with the CMB 
and Extra-galactic Background Light (EEL). In some models, effects of galactic and extra-galactic 
magnetic fields are considered. 

Dip Model The observed energy spectra of HiRes and TA above 10'** eV, in terms of the spectral 
shapes and their break point energies, are well explained by the GZK process; interactions of extra- 
galactic protons with the CMB producing the "ankle" via e^ e" pair production at around 10'^^ eV, 
and the "cutoff" via pion photo-production around 10'^^ eV (see Fig. |2]by R. Aloisio). The energy 
Ei/2 of the integrated flux is lo''^^^*""^ eV for the HiRes spectrum, and agrees with the theoretical 
calculation of lO'**'^^ eV {V. Berezinsky) within a systematic uncertainty of the energy scale of HiRes 
(17%). This picture (dip model) is consistent with the observation of a proton-like composition above 
10'** eV by HiRes and TA using the ji,„ax measurement. 

Multi-Composition Model On the other hand. Auger, with its measured composition changing 
from proton to heavier nuclei in this energy range, has not associated the observed spectral features 
with a particular physics mechanism. A variety of models with different source composition, accelera- 
tion and propagation assumptions have been proposed to fit Auger's high statistics spectrum data, and 
measurements of < X„u,x > and RMS(X„,fl^). See Fig. [3] for example for the case of p-He-Fe multi- 
composition by R. Aloisio. Although these models often contain CR acceleration hmits at the origin as 
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important ingredients, it is generally conceived that the photo-disintegration of iron nuclei by the CMB 
is the dominant cause creating the high energy flux suppression. The effect of photo-disintegration on 
nuclei was already mentioned in the original GZK papers. The ankle may be explained in these models 
as the acceleration limit of extra-galactic protons, or of the galactic heavy components corresponding 
to a galactic to extra-galactic transition of CRs. 





E(eV) 

Fig. 2. HiRes spectrum and the prediction of the dip Fig. 3. Auger spectrum and prediction by one of the 
model (R. Aloisio) multi-composition models {R. Aloisio). 



The following experimental challenges are identified to pin down the correct physics model of 
UHECR spectrum: 

- Energy scale uncertainty in FD: The flux suppression by the GZK process, whether it is for 
proton or for iron, is absolute in energy scale, and the parameter such as E 1/2 is relatively insensitive 
to the choice of the spectral index at the source and the cosmological evolution of CRs. Note that, 
between Auger and TA, there is a significant difference in the energy scale; Auger is ~18% lower 
than TA in the rescaling of energy spectra, and Auger's cutoff breakpoint energy is ~46% lower 
than that of TA (lO'''-^^ eV for TA and lO''' '*! eV for Auger See [9]). 

- Energy linearity in SD: The calibration of SD energy is made only by the limited sample of 
SD/FD coincident events, i.e. the direct calibration is limited to energies below the onset of the 
high-energy cutoff. The extrapolation to higher energies relies on the validity of the Constant In- 
tensity Cut (CIC) method for Auger, and Monte Carlo air shower simulation for TA, both of which 
require independent experimental verification. 

- Mass Composition: The mass composition plays a key role in distinguishing different models. 
The change of < Xj^ix > also directly affects the energy determination. Meaningful measurement 
of Xmo j. is now limited to 10"^" eV - 10'^'^ eV. An extension of composition measurement to both 
higher and lower energy regions is important. 

- Around the cutoff and above: The spectral shape above 10'^^ eV may be a convolution of the 
GZK cutoff and the acceleration limit of UHECRs. The cutoff shape also carries information on 
the composition as seen in Figs|2]and|3] Above 10'"^^ eV, the agreement of spectral shapes among 
experiments is not as satisfactory as the lower energy region (see FiglTJ. In order to determine the 
exact shape of the cutoff above 10'^'^ eV, measurements with higher statistics and less systematic 
error are essential. 



6 Mass Composition 

The primary mass composition of UHECRs is determined from the X,„ax measurement. The X„iax can 
be determined for each FD event, but the mass composition is usually determined only statistically 
by comparing the < X,„a_^ > and RMS(X,„flv) with the expectation from the EAS simulation. This 
is due to the large event-by-event fluctuation of X,„ax, particularly for proton events, and the limited 
accuracy of X„,ax measurement. The most straight-forward way of measuring X„,ax is by imaging the 
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EAS development in the atmosphere using the FD, which has been adopted by HiRes, Auger and TA. 
There are less direct way of measuring the Cherenkov light distribution on the ground (by Yakutsk) 
and muon detection on the ground (by Auger/SD). 

Measurements of < X„a.i > and RMS(X,„aj) are compiled by the composition WG 1 10|. They are 
shown in Fig. |4] for Auger and Yakutsk, in Figs. |5] and |6] for HiRes, in Fig. Q for TA (preliminary) 
together with the results of air shower simulations for protons and iron. The different measurements 
are not overlaid in one plot because each experiment treats, either includes or removes, analysis biases 
and resolution effects differently, and no direct comparison of measured values is adequate. 




Fig. 4. The energy dependence of < ^,„nv > for Auger and Yakutsk (left). The measured < X„„_^ > for Auger 
and Yakutsk has negligible analysis biases, and is directly comparable with the simulation (red and blue lines). 
The energy dependence of RMS(X,„av) for Auger and Yakutsk is shown on the right. The RMS(X,„av) of Auger 
is defined as the r.m.s. of the X„,nj- deviations. The experimental resolution estimated from the simulation is 
subtracted in quadrature. The highest energy bin of Yakutsk contains only 3 events. 
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Fig. 5. The energy dependence of < X,„„y > for HiRes. 
The measured < Z^a,- > is compared with the simula- 
tion (lines) with analysis biases included both for the 
measurement and simulation. 
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Fig. 6. The energy dependence of RMS(X„,„j-) for 
HiRes. The RMS(Xmaj) is defined as the cr obtained 
by fitting the X„nj. distribution with a Gaussian trun- 
cated at 2(T. The experimental resolution is included 
for both the measurement and simulation. 



As seen in the plots, all experiments agree that the light component (proton) dominates in the en- 
ergy region below the ankle for 10'**'^ eV - 10'^'^ eV. Above 10"*'^ eV, Auger's measurement suggests 
a change of composition from light to heavier, reaching almost iron at 10'^^ eV. This is seen in both 
< Xmax > and RMS(X,„flf), but is more pronounced in RMS(Xmfl.v)- On the other hand, no clear indica- 
tion of the composition change is seen for the < X„„„ > of HiRes, TA and Yakutsk up approximately 
to lO'''-^ eV. 

Distinctions between Auger and HiRes/TA/ Auger is obviously that Auger covers the southern sky, 
and HiRes/TA/Yakutsk cover the northern sky. Among other important distinctions are event detection 
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Fig. 7. The energy dependence of < X„,„y > for TA 
(preliminary). The measured < X„,a.x > is compared 
with the simulation (lines) with analysis biases in- 
cluded for both the measurement and simulation. 
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Fig. 8. The energy dependence of average composition 
<lnA> using QGSJET-11. The shaded area indicates sys- 
tematic uncertainties of each experiment. The energies 
are rescaled as in the case for the spectrum comparison. 



topologies and data analysis methods; Auger uses hybrid events using FD and SD (muon based) in 
coincidence, and HiRes and TA use stereo FD events using two FDs (EM energy based) in coincidence. 
In the data analysis. Auger uses tight acceptance cuts to minimize the analysis bias, whereas HiRes 
and TA use looser acceptance cuts to accept more events and use elaborate simulations to estimate the 
acceptance bias. The composition WG scrutinized the analysis methods of Auger, HiRes and TA in 
order to understand whether the analysis method is responsible for the difference. 

In order to guarantee good measurements of X„,a j, all experiments apply the so-called bracket cut 
for the acceptance. This cut requires that the of individual events be observed in the Field of 
View (FoV) of the FD telescope. The tight bracket condition, however, causes certain bias in < X,nax > 
and RMS(X,„aj:) because, for a given event geometry, we do not know a priori how many unbracketed 
events exist outside the FoV. 

HiRes and TA developed an elaborate EAS event simulation, analyzed the generated MC events 
using the same calibration and analysis programs as the data, and the data-MC agreement was checked. 
The event reconstruction was made by the inverse-MC method using the same parameter set for the 
data and MC. The data observables, including < X,„„i > and RMS(X,„fl^) with its energy dependence, 
are well reproduced by the pure proton primary using the QGSJET-II interaction model, and this is the 
basis for the claim of an unchanged proton component by HiRes and TA. The estimated acceptance 
bias of TA is less than 15 g/cm^ for QGSJET-II proton above 10"^ eV, and is negligible for QGSJET- 
II iron. The reconstruction bias, defined as the difference between the real (thrown) X,„fl^ and the 
reconstructed X„,af, is approximately 10 g/cm^ independent of the energy and the composition. 

Auger's X„,a , analysis requires tighter condition for the acceptance such that an event is accepted 
only if not only the observed X,„a f but also an entire range of the probable X„,a j of that event geometry 
(and energy) is included in the FoV. The range of accepted Y^max is determined by checking that the 
resultant < X„,a j > is not influenced by the FoV cuts IJOJ . This method is tuned to minimize the Y^max 
acceptance bias, at the cost of fewer accepted events; the tight FoV cut reduces the accepted events of 
Auger by half for energies above 10'^^ eV. The performance of the analysis was estimated by using 
the MC; the combined (reconstruction + acceptance) X,„fl^ bias is smaller than 4 g/cm^, and the X„,fl j 
resolution is 20-25 g/cm^. 

The composition WG made numerous checks on the analysis results, and the data-MC agreement 
was examined for many observables as reported in [10|. In shoit, no particular flaw of analysis, or 
definite clue for the resolution of the difference was identified so far. The WG is now proceeding to 
the next step to analyze the Auger data with TA's looser cuts, to analyze the TA data with Auger's tight 
cut, in order to see whether the result is unchanged. Work is also in progress to generate and analyze 
MC event sets having the same y^max distribution as observed by Auger, and pass them through the 
analysis chain of TA and HiRes. 
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The results of < X„iax > from different experiments (or analysis methods) may be compared by 
converting the < X„,ax > into <lnA> flOl, as it is known that the change of X,„„x is approximately 
proportional to InA at a fixed CR energy, where A is the atomic number of the primary CR. The 
composition WG calculated the value of <lnA> for each experiment by proportionating the difference 
of the measured X,„af between the expected X„„^ for protons and iron. The result obtained by using 
QGSJET-II is shown in Fig. [8] The shaded area indicates the systematic errors of each experiment; 
12 g/cm^ for Auger and TA, 20 g/cm^ for Yakutsk, 6 g/cm^ for HiRes. It is seen that HiRes and 
Auger give different results and are not compatible for energies above 10 eV. The results of TA and 
Yakutsk partially overlap both with HiRes and Auger A comparison using the SIB YLL model makes 
the difference between Auger and HiRes smaller, but they do not overlap above 10'**^ eV. 

It should be noted that the statistics of Auger is higher than all other experiments combined, even 
with the tight acceptance cutQ. The WG report demonstrated that the energy dependence of Auger's 
<lnA> is statistically not consistent with a constant value of <lnA>, whereas the same plots by other 
experiments are statistically consistent with both constant and changing values of <lnA> with respect 
to log(E) (see Fig. 12 of [ 10|). More events (preferably x3 or more) in future, particularly in the north- 
ern hemisphere, are needed to conclude definitively whether the composition stays light as the energy 
increases, or turns heavy as in the case of Auger in the southern hemisphere. In addition, on-going 
improvements of HAS simulation with better models and newly available collider data (see section |9|l 
are expected to contribute a lot to predict the < X„ax >, RMS(X„„„) and its energy development in a 
more reliable way. The composition WG will continue its activity in future. 

Finally, the FD and SD hardware and their calibration needs due attention for the X„tax comparison. 
This becomes important when the comparison of analyses does not solve the problem. HiRes, TA 
and Auger use different optics, PMT matrix and signal recording methods for the FD telescope. The 
atmospheric models and its attenuation effects are different for each experimental site. S. Ogio, in his 
talk of future plans Il20l . proposed to exchange a set of SDs and FDs between Auger and TA sites, 
and to make a simultaneous observation of the same UHECR event. This will be very effective in 
excavating systematics in the detector and analysis, if any important one exists. 



7 Anisotropy 

It is expected that the magnetic deflections of UHECRs by the galactic and extra-galactic magnetic 
fields decrease down to a few degrees for energies larger than 10^" eV, and a correlation between 
the UHECR arrival direction and the source "star" would be showing up. This correlation would be 
particularly distinct if the source astrophysical objects are limited within the GZK horizon (~100 
Mpc). In 2007, Auger published a correlation of UHECR (E > 57 EeV|f|) arrival directions and nearby 
AGNs (z < 0.018, or within 75 Mpc) in the Veron-Cetty and Veron (VCV) catalog, 12th edition within 
a separation of 3.1° [14|. The result was updated in the Beijing ICRC conference and its skymap 
and time development diagram are shown in Figs. |9]and[T0] The rate of correlated UHECR events 
presented in Beijing decreased from 61% to 33% but it is still ~3.0cr (P=0.006) away from a chance 
correlation. A concentration of correlated events in the vicinity of Cen A was reported in I.15J , but the 
significance did not increase despite higher statistics lfT6l . 

HiRes searched for a correlation in the northern hemisphere but did not find any. TA searched for 
correlations in the recent SD data sample with zenith angles less than 45°, and found 11 events are 
correlated from the total of 25 events under the same conditions as the Auger AGN correlation ifTTI 
0. The number of correlation expected from the uniform and random set of 25 events is 5.9, or the 
observed "correlation" is more than 2cr (P=0.02) away from the accidental correlation. Both HiRes 
and TA did not optimize the search parameters and used the same conditions as Auger in order to stay 

'* The number of events above 10" eV is 452 for Auger, 123 for HiRes, 67 for TA and 22 for Yakutsk. 

^ 57 EeV = 10" '* eV, EeV = 10'^ eV. The original publication had a threshold energy of 56 EeV. It was later 
stated to be equivalent to 57 EeV in the updated energy scale. 

* The energy scales of TA/HiRes may be ~20% larger than that of Auger as seen in section|4] and the effective 
energy threshold for the search may be different by this amount. 
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free from the search bias. The skymap of TA and the time development of the correlation are shown in 
Figs.[n]and[l2] 
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Fig. 9. Arrival directions of UHECRs (E > 57 EeV) 
in tlie souttiem sky observed by Auger. A correlation 
with nearby AGNs (z < 0.018) was suggested. 
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Fig. 10. Time development of Auger AGN correlation. 
The ratio of AGN-correlated events over the total is 
plotted. The dashed line corresponds to the ratio ex- 
pected from the accidental coincidence. 
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Fig. 11. Arrival directions of UHECRs in the northern 
sky observed by TA. The same correlation conditions 
as in Fig.|9]are applied. 
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Fig. 12. Time development of AGN correlation of TA. 
The number of AGN-correlated events is plotted. The 
solid line corresponds to the number expected from the 
accidental coincidence. 



Auto-correlations among UHECR events are searched for by HiRes, Auger and TA but no positive 
results have been reported. A global harmonic analysis was performed by Auger, and correlations with 
the Large Scale Structure (LSS) of galaxies have been searched for by HiRes, Auger and TA, but the 
results are so far consistent with both the isotropy and the LSS. 

The anisotropy WG summarized that, despite a large effort spent on the search for anisotropy and 
correlation, no clear signal has so far been established with certainty [ 1 1 1. Various hints, however, are 
emerging from the observations of Auger and TA, and they should be closely followed in coming days. 
The prospects and challenges by the anisotropy WG are reiterated as follows: 

- Composition: Expectation for correlation obviously depends on the actual composition, or the 
atomic number, of the primary UHECRs. It should be noted that above 10'^^ eV where the corre- 
lations with AGNs and LSS are expected, no direct measurement of composition by X„a.i exists. 
Discrimination power for composition on an event-by-event basis would be a great asset for the 
future anisotropy search. 

- Large Exposure: The next 5 years will be an exciting time when some of the hints of anisotropy 
seen by Auger and TA will be claiified with more statistics. In order to fully test the observed 
anisotropy, or to find new kind of anisotropy and characterize its cause, one order of magnitude 
larger exposure is needed by the next generation observatory. 
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- Energy Assignment: An accurate energy assignment of each event is essential for making energy 
cuts for trans-GZK anisotropy and also for large-scale anisotropy studies at EeV energies. The 
energy scale error (of FD) is expected to be improved from the present ~20% down to ~10%, 
matching the energy resolution (of SD) at the highest energy region. The present angular resolution 
of ~1° may be sufficient for the future anisotropy study Q 

- Full-sky coverage: There are many advantages to operate an observatory with full-sky coverage; 
all possible point sources would be exposed, and the structure of anisotropy with minimal expo- 
sure distortion could be achieved. Besides, with full-sky coverage, the methods of spherical har- 
monics and multipole analysis could be applied without additional assumptions on the anisotropy 
behavior in the uncovered region of the sky. Unified analysis of southern and northern skies are 
common interests of existing groups; such analysis was initiated using the Auger and TA data by 
the anisotropy WG, and is being pursued now. 

- Wide Energy Coverage: Currently preferred models of UHECR (see section |5) suggest that 
cosmic rays are extra-galactic above 10'^ eV, but the scenario of transition from galactic to extra- 
galactic below 10'** eV is not established, and the predictions of models differ (see Berezinsky). 
Understanding this transition is an important step forward to provide a complete picture for origins 
of UHECRs. The transition should be signaled by the composition change (from heavy to light), 
spectrum deformation (second knee?) and the energy-dependent anisotropy as hinted by the phase 
change in the harmonic analysis IfTTI of Auger Such low-energy extensions are proceeding well 
as HEAT and AMIGA in Auger and TALE in TA to supply high quality, high statistics data down 
to 10'^'^ eV in the near future. 



8 Multi-Messenger (UHE gammas and neutrinos) 

Any top-down model of UHECR origin expects abundant UHE y rays and neutrinos as a decay product 
of n mesons. The GZK process also produces UHE gammas and neutrinos from pion decays of the 
excited nucleon resonance. No such UHE cosmic y's and y's have been observed so far by the UHECR 
experiments. 

Limits on the flux of UHE y's by Auger, TA and Yakutsk already place strong constraints on the 
top-down models as seen in Fig. [T3]lfT2l . The multi-messenger WG further compiled the present (and 
"anticipated" by 2015) experimental upper limits and compared them with the expected flux from 
cosmogenic (GZK) production models. They are reproduced here in FigsfT4land[T5] 




E„l„. eV 

Fig. 14. Experimental limits placed on the UHE pho- 
Fig. 13. Experimental limits placed on the UHE pho- ton flux. Expected limits by 2015 are indicated by 
ton flux and expectation from several top-down accel- lines. The shaded region corresponds to the flux pre- 
eration models of UHECRs. dieted by several cosmogenic models 11121 . 



^ For TA/SD events, the energy resolution is 13% and the angular resolution is 1.1° for E > 57 EeV by simula- 
tion studies. Other experiments have similar numbers. 
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Fig. 15. Experimental limits placed on the UHE neutrino flux. Expected fluxes from several cosmogenic models 
are shown by the shaded region and black lines. For details, see WG's report Iil2il . 



The expected sensitivities for UHE y's in Fig.[T4]tells us that the next generation UHECR detector 
with one order of magnitude larger exposure than the present Auger and TA has a significant chance to 
detect cosmogenic y rays, provided the primary composition is proton. The iron primary is expected 
to give an order of magnitude smaller flux than the proton primary. The UHECR detectors planned 
for future, such as the nGO and JEM/EUSO (section [TOli may be the only detectors to observe UHE 
cosmogenic y rays. 

For the case of cosmogenic v's, as seen in Fig [15] the recently completed IceCube (IC-86) will 
have a good chance of detecting UHE v's after ~5 years of running in the region of 10^^ eV - 10'^ 
eV, for the case of primary protons. The south pole radio experiments, ANITA by balloon and RICE 
in the ice, already contribute to set upper limits for energies greater than 10'^ eV. Future huge ground 
array nGO with 10 times larger acceptance than the present Auger, covers the UHE v detection in the 
region of 10'^ eV - 10^° eV, and in case if detected, it would give most precise information on the 
energy and arrival direction, and possibly the species (e, //, r) of the UHE neutrinos. The JEM/EUSO 
will have the largest sensitivity but its energy acceptance is limited above 10'''^ eV. The iron primary 
or the mixed composition may reduce the detection possibility of cosmogenic neutrinos by roughly an 
order of magnitude. 



9 Hadronic Interaction and AS simulation 

Despite the success of determining the UHECR energy by the new air fluorescence technique, and 
identifying the GZK traces in the energy spectrum, understanding the EAS itself is still incomplete 
and unsatisfactory. The problem is demonstrated in the following discrepancies between the EAS 
observation and the predictions of EAS simulation: 

- Number of Muons: The number of muons measured on the ground, which is a tracer of hadronic 
energies in the EAS, is up to ~100% larger than that predicted by the standard air shower simula- 
tion (see A. Yushkov and G. Rodriguez). It is certain that hadronic interaction models at UHE, as 
well as the nuclear elfect, is responsible for this. Incomplete hadronic interaction database at lower 
energies may also be contributing to this discrepancy. 

- Longitudinal Shower Development: Different interaction models give different longitudinal de- 
velopment and X,„fl V distribution. The prediction of < X,„ax >, RMS(X,„fl v) and their energy devel- 
opment (elongation rate) differ among models. The measured < X^ax > and RMS(X,„a t) are not 
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well reproduced by the simulations. Cases of good agreement exist, such as the HiRes/TA mea- 
surements and QGSJET/proton simulation, but this does not sufficiently prove that "QGSJET" is 
the correct hadronic interaction model at UHE and the proton is the unique primary "composition" 
of UHECR. It is interesting to note, however, that existing model predictions for iron have little 
differences for < X,„ax >, RMS(X,„fl^) and its energy development. 

- Lateral Shower Development: The energy of UHECR events measured by the TA/SD is 27% 
higher than the energy measured by TA/FD. The AGASA/SD (plastic scintillator) gives higher 
energy than TA/SD, but only by about 9%. The systematic uncertainty of energy measurement is 
21% for TA/FD and -20% for TA/SD and AGASA, meaning that the difference is at the limit of 
both systematics. This difference means that the same EM energy measured 600~800 m away from 
the shower core by TA/SD and AGASA/SD is larger than the EM energy measured by the FD near 
the shower axis, suggesting the data tends to have larger lateral development than the simulation 
predicts. 

The advent of the LHC is opening up a possibility to improve this situation drastically. The LHC 
maximum beam energy is scheduled to be 7.0 TeV, corresponding to a lO'^ eV cosmic ray proton 
colliding with a proton target at rest. Hadronic and nuclear interactions at UHEs can be studied now 
at accelerators and the results can be fed back to the hadronic interaction models. Two very forward 
experiments at the LHC are particularly relevant. One is the LHCf measuring the spectrum of neutral 
particles (y, 7t°, and neutron) at 0°. The LHCf is dedicated for cosmic ray related studies. Another is 
the TOTEM experiment to measure total and (in-)elastic hadronic cross sections in the very forward 
region. It also offers valuable data to determine the development of EAS in the atmosphere. 

In addition to the heavy ion collision data available from the RHIC the LHC is foreseeing proton- 
lead colUsions next year, which will contribute to the understanding of the nuclear effect at UHE. This 
situation was reviewed and discussed by Y. Itow, T. Sako, T. Pierog, R. Engel and other symposium 
contributors. A report of hadronic interaction and simulation WG is also available 1 13 1. 

Using newly available data from accelerator experiments, an update of hadronic interaction models 
ai-e on the way for QGSJET (from v03 to v04) and EPOS (from vl.99 to -vLHC). T.Piemg reported 
the status of tuning with highlights on the following observables. 

- Total Cross Section: The total inelastic cross section determines the average depth of the first 
interaction of UHECRs in the atmosphere, and affects the speed of shower development and its 
fluctuation. The newly available data from TOTEM is being used for the tuning. 

- Multiplicity: This parameter, as well as the inelasticity, affects the speed of shower development, 
and the number of muons on the ground. The data from CMS and ALICE at large angles (small 
pseudo-rapidity) are being used for the tuning. 

- Inelasticity, or Energy Spectrum at Forward Region: The higher the leading particle energy, the 
slower the shower development expected. The data from CMS, ALICE and PHENIX are checked 
with the predictions of simulations. The p and pj distribution of y's and 7r°'s obtained by LHCf 
will be useful for this tuning. 

- (anti-)Nucleon Production at Very Forward Angles: The forward-produced nucleon is effective 
to maintain hadronic energies in the longitudinal shower development, thus affecting the muon rate 
at the ground level. The neutron production rate being analyzed by the LHCf experiment is eagerly 
waited. 

Auger recently published a measurement of the p-air total cross section using the X^ax tail distri- 
bution around 10'^ eV - 10'**^ eV where the dominant CR composition is considered to be proton 
|,18J . In future, with the progress of understanding the CR composition, a possible appearance of the 
new phase of hadronic interactions and/or unknown collective nuclear effects at UHE may be revealed 
in the data of EAS development. When it happens, it would become a good case for the bi-lateral 
contribution between the CR physics at UHE and elementary particle physics. 

It should be noted that the energy region above 10'^ eV is where we expect the galactic to extra- 
galactic transition of cosmic rays to take place (see section |7]i. On-going low energy extensions of 
HEAT/ AMIGA and TALE will be soon producing high statistics spectrum and composition data in 
this energy range, where the EAS simulation is being calibrated directly by the precision accelerator 
data. 
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10 Prospects 

A discussion on the future global facilities and the associated R&D were made in the symposium, 
keeping in mind the rapid progress of UHECR research and the still remaining experimental challenges 
as described above. 

- Next Generation Ground Observatory: Three conceptual designs for the next generation Ground 
Observatory with huge aperture (called nGO in this report) were presented by A. Letessier-Selvon, 
P. Privitera and 5'. Ogio. Astrophysical objectives of nGO were briefed by E. Parizot, and interdis- 
ciplinary science was discussed by L. Wiencke. 

- Space Observatory: Observation of UHECRs from space was reviewed by T. Ebisuzaki, and 
the mission, performance and the prototype test of JEM/EUSO were presented by A. Santangelo, 
M. Bertaina and M. Casolino. 

- R&D for Radio Detection: Radio detection of UHECRs at MHz range was discussed by O. Scholten 
and A. M. van den Berg. It was extensively studied by LOPES and CODALEMA, and now by 
AERA at the Auger site. C. Williams reported on the GHz radio measurement by the accelera- 
tor beam (MAYBE). Field tests for detecting the air shower radio signal at the Auger site (AM- 
BER, EASIER and MIDAS), and at KASCADE-Grande (CROME) were reported by P. Facal and 
R. Smida. J. Belz reported a plan of bistatic radar detection of UHECRs (TARA) at the TA site. 




Fig. 16. Footprint of the JEM-EUSO FoV pro- 
jected on the ground: nadir mode (green) and 
tilted mode (orange) are shown ||19] . 
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Fig. 17. The covered ground areas of past, present and future 
ground arrays. The proposed huge array by S.Ogio 11201 has 
10,000 TA-type plastic scintillator SDs deployed with 2.0 
km spacing and covers an area of 39.2k km^. 




Fig. 18. A possible configuration of the multi-component UHECR Observatory proposed by A. Letessier-Selvon 
1211 . It uses updated water tank SDs of Auger-type. Vertical segmentation of the tank serves for the EM calorime- 
try (upper thinner layer) and the muon counting (lower part). 
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A majority of the designs and R&D presented in the symposium concentrated on realizing a huge 
detection area, which is an order of magnitude larger than that of Auger. New generation experiments 
aim to collect ~1,000 events each year in the trans-GZK energy region (E > 10^'-^ eV) to identify the 
astrophysical origins of UHECRs. Their features may be compared as follows: 
Ground vs Space: The SDs of Auger and TA measure the UHECRs with an angular resolution of 
~1 .0° and an energy resolution of ~10% at 10^" eV. The nGO plans to extend the aperture maintaining 
the same level of resolution for trans-GZK events. This can be achieved by tuning the SD detector 
size and the separation using the presently available UHECR samples. A typical aperture of nGO is 
55,000 km^ sr (S. Ogio). For Auger and TA, the determination of composition by FD is statistical and 
is limited to hybrid events (duty < 10%). For nGO, a separate measurement of the EM energy and 
the number of muons particularly at each nGO/SD will give us a handle to determine the composition 
event-by-event (A. Letessier-Selvon). 

The JEM/EUSO aims at 2.5° and 30% resolution at 10^° eV. The '^,nax resolution is 120 g/cm^ and 
allows a separation of y's and v's from protons and iron. The energy threshold is ~10^'-^ eV. A large 
aperture of 100,000 km^ sr (after 20% duty factor applied) is obtained by a single telescope looking 
down on the earth's atmosphere from the International Space Station (ISS) at a height of 400 km. A tilt 
mode will increase the aperture to 350,000 km^ sr but with a sacrifice in performance. The JEM/EUSO 
design is tuned to search for astrophysical origins in the GZK horizon, and the all sky-coverage in the 
ISS orbit is advantageous for this purpose. 

Fluorescence vs Radio: A composition measurement by FD continues to be indispensable in the 
future nGO. The entire nGO area may be covered by an array of wide angle FD telescopes deployed in 
a mesh of ~20 km. Given that the event geometry is reconstructed using the nGO/SD information, only 
the time development of nGO/FD signal is sufficient to reconstruct the EAS longitudinal development. 
Such information can be collected by a small number of PMTs with much larger pixel size than that of 
Auger and TA, thus drastically reducing the cost for construction and maintenance {P. Privitera). The 
trigger capability, signal to noise ratio and resolution of '^max reconstruction are being studied for this 
scheme. 

When a radio detector is integrated in nGO, it would supply information of shower longitudinal 
development with high duty (~100%), httle atmospheric or weather correction and at low cost. Such 
a radio detector may be used in the nGO in place of the now standard FD telescopes. Many laboratory 
and field measurements have been performed to detect MHz, GHz and radio echo signals from EAS. 
The existence of MHz signal is firmly estabhshed by earlier experiments of LOPES, CODALEMA 
and others, and it is known that the signal emission is forward peaked and has a clear correlation with 
the geo-magnetic field. The GHz signal is expected to have an isotropic and un-polarized emission 
from the molecular bremsstrahlung, but so far the observed GHz signals by CROME and EASIER 
seem forward peaked, and not identified as the molecular bremsstrahlung signal. The forward peak, 
in both MHz and GHz cases, means the footprint of the signal on the ground is small and the radio 
measurement would be costly, requiring many antennae to cover the ground. The TARA radio echo 
test experiment is now set up on the TA site. It consists of a 54 MHz, 40 kW television transmitter 
illuminating the air above TA, and a set of radio receivers on the far side of TA, aiming to detect a 
forward scattered radio signal in coincidence with the TA measurement. 

So far, no proven radio method to replace the air fluorescence measurement is found, nevertheless, 
the R&D of radio signal detection continues in order to understand the nature of signal generation, and 
with its application for the X^a^ measurement in mind. 

Tank vs Scintillator: The Cherenkov signal of Auger water tanks is larger for penetrating muons 
than soft electrons and gammas in the shower, thus making the Auger/SD more sensitive to a change 
of composition with energy. Muons are also a good measure to investigate hadronic interactions in the 
shower development. 

The plastic scintillator of TA equally samples the muon and electron, and the resultant signal is 
dominated by the outnumbered electrons, or the EM energy. Since the bulk of the primary energy of 
UHECRs is transferred to the EM energy, separate measurement of the EM and muon components 
in the nGO/SD would allow us to take both advantages: a good linearity in the primary energy mea- 
surement by the EM component, and a good sensitivity to the primary composition (and hadronic 
interaction) by the muon counting. 
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For UHE 7 and v detection, a high profile of the water tank would allow a detection of very inclined 
EAS down to lower energies than thin plastic scintillators can achieve. The contrary may be true for 
the detection of UHE y ray interactions which contain much fewer muons than normal EAS of the 

same primary energy. 

In order to exchange ideas and experiences for the next generation detector, it was agreed to form 
a new working group for future projects. The JEM/EUSO group also expressed an interest to join the 
working group. The design of the next generation detector becomes more sensible and realistic when 
the current issues of primary composition and the energy scale/linearity are understood. An exchange 
of analysis method, calibration, simulation, obtained data and eventually the FD/SD detectors may be 
considered among experiments to attack this formidable problem. It is encouraging, as a first attempt 
of many others to be considered. Auger and TA initiated a photometric and optical calibration of FDs 
by flying a common standard light source in Utah and Malargue by using a GPS controlled octo-copter. 
The first flight at TA was completed in October 2012, and the next flight at Auger site is planned in 
November. 



1 1 Summary and Acknowledgement 

The scene of UHECR is rapidly evolving. We just settled a long lasting, difficult experimental question 
on the energy spectrum of extra galactic cosmic rays; Yes, there is a "cutoff" at 10^'-^ eV and "ankle" 
at 10'^^ eV. The accuracy of energy measurement is about 20% now by the new air fluorescence 
method, and it is steadily improving. In the northern hemisphere, a constant particle composition 
of proton is measured above 10'**'^ eV, whereas the composition is changing from proton to iron in 
the southern hemisphere. We eagerly examine our detectors and data analysis to know whether this 
difference is the nature or an experimental artifact. Calibration of air shower simulation is proceeding 
with the LHC and other accelerator data. We do hope it will shed light on the interpretation of X^ax 
data. No UHE gammas and neutrinos were detected, and the top-down acceleration of UHECR is 
strongly disfavored. Limits of astrophysical acceleration of UHECR is not seen, or may have been 
seen, depending on the composition result and its model interpretation. A transition of cosmic ray 
sources from galactic to extra-galactic may be identified in the near future by low energy extensions of 
present experiments (HEAT/ AMIGA and TALE). Signatures of UHECR sources come and go, though 
certain levels of anisotropy seem to be an unavoidable consequence of the present observations. By 
JEM/EUSO, space may become part of our work field in the near future. It will seek to identify 
UHECR origins in all sky in the trans-GZK region. The R&D and design of nGO, a new Ground 
Observatory, is proceeding. It will accumulate 1,000 events of well observed UHECRs each year to 
definitively solve existing questions on the UHECR, and to explore a new field of physics in the 
universe. So, be patient (P. Priviterd), and stay tuned. 

This summary report is based on the work of UHECR-2012 contributors and working groups. I 
have to regretfully admit the introduction of their works to this report is very incomplete. For more 
details and to understand complicated points, I encourage readers to consult the original papers in the 
same proceedings. To make this summary report, I am largely indebted for the work of contributors 
and WGs, but also feel responsible for my views expressed in this report, and any mistake, bias and 
inaccuracy appearing in this report. The UHECR-2012 was a very active and enjoyable symposium. 
High scientific activities and good spirits of the participants, and the symposium organizers, are truly 
appreciated. 
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